Genotypic and phenotypic analyses of a Pseudomonas aeruginosa chronic bronchiectasis isolate reveal differences from cystic fibrosis and laboratory strains

Background Pseudomonas aeruginosa is an environmentally ubiquitous Gram-negative bacterium and important opportunistic human pathogen, causing severe chronic respiratory infections in patients with underlying conditions such as cystic fibrosis (CF) or bronchiectasis. In order to identify mechanisms responsible for adaptation during bronchiectasis infections, a bronchiectasis isolate, PAHM4, was phenotypically and genotypically characterized. Results This strain displays phenotypes that have been associated with chronic respiratory infections in CF including alginate over-production, rough lipopolysaccharide, quorum-sensing deficiency, loss of motility, decreased protease secretion, and hypermutation. Hypermutation is a key adaptation of this bacterium during the course of chronic respiratory infections and analysis indicates that PAHM4 encodes a mutated mutS gene responsible for a ~1,000-fold increase in mutation rate compared to wild-type laboratory strain P. aeruginosa PAO1. Antibiotic resistance profiles and sequence data indicate that this strain acquired numerous mutations associated with increased resistance levels to β-lactams, aminoglycosides, and fluoroquinolones when compared to PAO1. Sequencing of PAHM4 revealed a 6.38 Mbp genome, 5.9 % of which were unrecognized in previously reported P. aeruginosa genome sequences. Transcriptome analysis suggests a general down-regulation of virulence factors, while metabolism of amino acids and lipids is up-regulated when compared to PAO1 and metabolic modeling identified further potential differences between PAO1 and PAHM4. Conclusions This work provides insights into the potential differential adaptation of this bacterium to the lung of patients with bronchiectasis compared to other clinical settings such as cystic fibrosis, findings that should aid the development of disease-appropriate treatment strategies for P. aeruginosa infections

Risks, benefits, and knowledge gaps of non-native tree species in Europe

Changing ecosystem conditions and diverse socio-economical events have contributed to an ingrained presence of non-native tree species (NNTs) in the natural and cultural European landscapes. Recent research endeavors have focused on different aspects of NNTs such as legislation, benefits, and risks for forestry, emphasizing that large knowledge gaps remain. As an attempt to fulfill part of these gaps, within the PEN-CAFoRR COST Action (CA19128) network, we established an open-access questionnaire that allows both academic experts and practitioners to provide information regarding NNTs from 20 European countries. Then, we integrated the data originating from the questionnaire, related to the country-based assessment of both peer-reviewed and grey literature, with information from available datasets (EUFORGEN and EU-Forest), which gave the main structure to the study and led to a mixed approach review. Finally, our study provided important insights into the current state of knowledge regarding NNTs. In particular, we highlighted NNTs that have shown to be less commonly addressed in research, raising caution about those characterized by an invasive behavior and used for specific purposes (e.g., wood production, soil recultivation, afforestation, and reforestation). NNTs were especially explored in the context of resilient and adaptive forest management. Moreover, we emphasized the assisted and natural northward migration of NNTs as another underscored pressing issue, which needs to be addressed by joint efforts, especially in the context of the hybridization potential. This study represents an additional effort toward the knowledge enhancement of the NNTs situation in Europe, aiming for a continuously active common source deriving from interprofessional collaboration. Copyright © 2022 Dimitrova, Csilléry, Klisz, Lévesque, Heinrichs, Cailleret, Andivia, Madsen, Böhenius, Cvjetkovic, De Cuyper, de Dato, Ferus, Heinze, Ivetić, Köbölkuti, Lazarević, Lazdina, Maaten, Makovskis, Milovanović, Monteiro, Nonić, Place, Puchalka and Montagnoli

Improving the iMM904 S. cerevisiae metabolic model using essentiality and synthetic lethality data

<p>Abstract</p> <p>Background</p> <p><it>Saccharomyces cerevisiae </it>is the first eukaryotic organism for which a multi-compartment genome-scale metabolic model was constructed. Since then a sequence of improved metabolic reconstructions for yeast has been introduced. These metabolic models have been extensively used to elucidate the organizational principles of yeast metabolism and drive yeast strain engineering strategies for targeted overproductions. They have also served as a starting point and a benchmark for the reconstruction of genome-scale metabolic models for other eukaryotic organisms. In spite of the successive improvements in the details of the described metabolic processes, even the recent yeast model (i.e., <it>i</it>MM904) remains significantly less predictive than the latest <it>E. coli </it>model (i.e., <it>i</it>AF1260). This is manifested by its significantly lower specificity in predicting the outcome of grow/no grow experiments in comparison to the <it>E. coli </it>model.</p> <p>Results</p> <p>In this paper we make use of the automated GrowMatch procedure for restoring consistency with single gene deletion experiments in yeast and extend the procedure to make use of synthetic lethality data using the genome-scale model <it>i</it>MM904 as a basis. We identified and vetted using literature sources 120 distinct model modifications including various regulatory constraints for minimal and YP media. The incorporation of the suggested modifications led to a substantial increase in the fraction of correctly predicted lethal knockouts (i.e., specificity) from 38.84% (87 out of 224) to 53.57% (120 out of 224) for the minimal medium and from 24.73% (45 out of 182) to 40.11% (73 out of 182) for the YP medium. Synthetic lethality predictions improved from 12.03% (16 out of 133) to 23.31% (31 out of 133) for the minimal medium and from 6.96% (8 out of 115) to 13.04% (15 out of 115) for the YP medium.</p> <p>Conclusions</p> <p>Overall, this study provides a roadmap for the computationally driven correction of multi-compartment genome-scale metabolic models and demonstrates the value of synthetic lethals as curation agents.</p

A Method for Efficient Calculation of Diffusion and Reactions of Lipophilic Compounds in Complex Cell Geometry

A general description of effects of toxic compounds in mammalian cells is facing several problems. Firstly, most toxic compounds are hydrophobic and partition phenomena strongly influence their behaviour. Secondly, cells display considerable heterogeneity regarding the presence, activity and distribution of enzymes participating in the metabolism of foreign compounds i.e. bioactivation/biotransformation. Thirdly, cellular architecture varies greatly. Taken together, complexity at several levels has to be addressed to arrive at efficient in silico modelling based on physicochemical properties, metabolic preferences and cell characteristics. In order to understand the cellular behaviour of toxic foreign compounds we have developed a mathematical model that addresses these issues. In order to make the system numerically treatable, methods motivated by homogenization techniques have been applied. These tools reduce the complexity of mathematical models of cell dynamics considerably thus allowing to solve efficiently the partial differential equations in the model numerically on a personal computer. Compared to a compartment model with well-stirred compartments, our model affords a more realistic representation. Numerical results concerning metabolism and chemical solvolysis of a polycyclic aromatic hydrocarbon carcinogen show good agreement with results from measurements in V79 cell culture. The model can easily be extended and refined to include more reactants, and/or more complex reaction chains, enzyme distribution etc, and is therefore suitable for modelling cellular metabolism involving membrane partitioning also at higher levels of complexity

Phenotypic and Transcriptomic Response of Auxotrophic Mycobacterium avium Subsp. paratuberculosis leuD Mutant under Environmental Stress

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of severe gastroenteritis in cattle. To gain a better understanding of MAP virulence, we investigated the role of leuD gene in MAP metabolism and stress response. For this, we have constructed an auxotrophic strain of MAP by deleting the leuD gene using allelic exchange. The wildtype and mutant strains were then compared for metabolic phenotypic changes using Biolog phenotype microarrays. The responses of both strains to physiologically relevant stress conditions were assessed using DNA microarrays. Transcriptomic data was then analyzed in the context of cellular metabolic pathways and gene networks. Our results showed that deletion of leuD gene has a global effect on both MAP phenotypic and transcriptome response. At the metabolic level, the mutant strain lost the ability to utilize most of the carbon, nitrogen, sulphur, phosphorus and nutrient supplements as energy source. At the transcriptome level, more than 100 genes were differentially expressed in each of the stress condition tested. Systems level network analysis revealed that the differentially expressed genes were distributed throughout the gene network, thus explaining the global impact of leuD deletion in metabolic phenotype. Further, we find that leuD deletion impacted metabolic pathways associated with fatty acids. We verified this by experimentally estimating the total fatty acid content of both mutant and wildtype. The mutant strain had 30% less fatty acid content when compared to wildtype, thus supporting the results from transcriptional and computational analyses. Our results therefore reveal the intricate connection between the metabolism and virulence in MAP

Metabolic Networks of Sodalis glossinidius: A Systems Biology Approach to Reductive Evolution

Background: Genome reduction is a common evolutionary process affecting bacterial lineages that establish symbiotic or pathogenic associations with eukaryotic hosts. Such associations yield highly reduced genomes with greatly streamlined metabolic abilities shaped by the type of ecological association with the host. Sodalis glossinidius, the secondary endosymbiont of tsetse flies, represents one of the few complete genomes available of a bacterium at the initial stages of this process. In the present study, genome reduction is studied from a systems biology perspective through the reconstruction and functional analysis of genome-scale metabolic networks of S. glossinidius. Results: The functional profile of ancestral and extant metabolic networks sheds light on the evolutionary events underlying transition to a host-dependent lifestyle. Meanwhile, reductive evolution simulations on the extant metabolic network can predict possible future evolution of S. glossinidius in the context of genome reduction. Finally, knockout simulations in different metabolic systems reveal a gradual decrease in network robustness to different mutational events for bacterial endosymbionts at different stages of the symbiotic association. Conclusions: Stoichiometric analysis reveals few gene inactivation events whose effects on the functionality of S. glossinidius metabolic systems are drastic enough to account for the ecological transition from a free-living to hostdependent lifestyle. The decrease in network robustness across different metabolic systems may be associated with th

Genome-Scale Reconstruction of Escherichia coli's Transcriptional and Translational Machinery: A Knowledge Base, Its Mathematical Formulation, and Its Functional Characterization

Metabolic network reconstructions represent valuable scaffolds for ‘-omics’ data integration and are used to computationally interrogate network properties. However, they do not explicitly account for the synthesis of macromolecules (i.e., proteins and RNA). Here, we present the first genome-scale, fine-grained reconstruction of Escherichia coli's transcriptional and translational machinery, which produces 423 functional gene products in a sequence-specific manner and accounts for all necessary chemical transformations. Legacy data from over 500 publications and three databases were reviewed, and many pathways were considered, including stable RNA maturation and modification, protein complex formation, and iron–sulfur cluster biogenesis. This reconstruction represents the most comprehensive knowledge base for these important cellular functions in E. coli and is unique in its scope. Furthermore, it was converted into a mathematical model and used to: (1) quantitatively integrate gene expression data as reaction constraints and (2) compute functional network states, which were compared to reported experimental data. For example, the model predicted accurately the ribosome production, without any parameterization. Also, in silico rRNA operon deletion suggested that a high RNA polymerase density on the remaining rRNA operons is needed to reproduce the reported experimental ribosome numbers. Moreover, functional protein modules were determined, and many were found to contain gene products from multiple subsystems, highlighting the functional interaction of these proteins. This genome-scale reconstruction of E. coli's transcriptional and translational machinery presents a milestone in systems biology because it will enable quantitative integration of ‘-omics’ datasets and thus the study of the mechanistic principles underlying the genotype–phenotype relationship

Global maps of soil temperature

Research in global change ecology relies heavily on global climatic grids derived from estimates of air temperature in open areas at around 2&nbsp;m above the ground. These climatic grids do not reflect conditions below vegetation canopies and near the ground surface, where critical ecosystem functions occur and most terrestrial species reside. Here, we provide global maps of soil temperature and bioclimatic variables at a 1-km2 resolution for 0\u20135 and 5\u201315&nbsp;cm soil depth. These maps were created by calculating the difference (i.e. offset) between in situ soil temperature measurements, based on time series from over 1200 1-km2 pixels (summarized from 8519 unique temperature sensors) across all the world's major terrestrial biomes, and coarse-grained air temperature estimates from ERA5-Land (an atmospheric reanalysis by the European Centre for Medium-Range Weather Forecasts). We show that mean annual soil temperature differs markedly from the corresponding gridded air temperature, by up to 10\ub0C (mean&nbsp;=&nbsp;3.0&nbsp;\ub1&nbsp;2.1\ub0C), with substantial variation across biomes and seasons. Over the year, soils in cold and/or dry biomes are substantially warmer (+3.6&nbsp;\ub1&nbsp;2.3\ub0C) than gridded air temperature, whereas soils in warm and humid environments are on average slightly cooler ( 120.7&nbsp;\ub1&nbsp;2.3\ub0C). The observed substantial and biome-specific offsets emphasize that the projected impacts of climate and climate change on near-surface biodiversity and ecosystem functioning are inaccurately assessed when air rather than soil temperature is used, especially in cold environments. The global soil-related bioclimatic variables provided here are an important step forward for any application in ecology and related disciplines. Nevertheless, we highlight the need to fill remaining geographic gaps by collecting more in situ measurements of microclimate conditions to further enhance the spatiotemporal resolution of global soil temperature products for ecological applications

Global maps of soil temperature.

Research in global change ecology relies heavily on global climatic grids derived from estimates of air temperature in open areas at around 2 m above the ground. These climatic grids do not reflect conditions below vegetation canopies and near the ground surface, where critical ecosystem functions occur and most terrestrial species reside. Here, we provide global maps of soil temperature and bioclimatic variables at a 1-km &lt;sup&gt;2&lt;/sup&gt; resolution for 0-5 and 5-15 cm soil depth. These maps were created by calculating the difference (i.e. offset) between in situ soil temperature measurements, based on time series from over 1200 1-km &lt;sup&gt;2&lt;/sup&gt; pixels (summarized from 8519 unique temperature sensors) across all the world's major terrestrial biomes, and coarse-grained air temperature estimates from ERA5-Land (an atmospheric reanalysis by the European Centre for Medium-Range Weather Forecasts). We show that mean annual soil temperature differs markedly from the corresponding gridded air temperature, by up to 10°C (mean = 3.0 ± 2.1°C), with substantial variation across biomes and seasons. Over the year, soils in cold and/or dry biomes are substantially warmer (+3.6 ± 2.3°C) than gridded air temperature, whereas soils in warm and humid environments are on average slightly cooler (-0.7 ± 2.3°C). The observed substantial and biome-specific offsets emphasize that the projected impacts of climate and climate change on near-surface biodiversity and ecosystem functioning are inaccurately assessed when air rather than soil temperature is used, especially in cold environments. The global soil-related bioclimatic variables provided here are an important step forward for any application in ecology and related disciplines. Nevertheless, we highlight the need to fill remaining geographic gaps by collecting more in situ measurements of microclimate conditions to further enhance the spatiotemporal resolution of global soil temperature products for ecological applications

Global maps of soil temperature

Research in global change ecology relies heavily on global climatic grids derived from estimates of air temperature in open areas at around 2 m above the ground. These climatic grids do not reflect conditions below vegetation canopies and near the ground surface, where critical ecosystem functions occur and most terrestrial species reside. Here, we provide global maps of soil temperature and bioclimatic variables at a 1-km2 resolution for 0–5 and 5–15 cm soil depth. These maps were created by calculating the difference (i.e. offset) between in situ soil temperature measurements, based on time series from over 1200 1-km2 pixels (summarized from 8519 unique temperature sensors) across all the world\u27s major terrestrial biomes, and coarse-grained air temperature estimates from ERA5-Land (an atmospheric reanalysis by the European Centre for Medium-Range Weather Forecasts). We show that mean annual soil temperature differs markedly from the corresponding gridded air temperature, by up to 10°C (mean = 3.0 ± 2.1°C), with substantial variation across biomes and seasons. Over the year, soils in cold and/or dry biomes are substantially warmer (+3.6 ± 2.3°C) than gridded air temperature, whereas soils in warm and humid environments are on average slightly cooler (−0.7 ± 2.3°C). The observed substantial and biome-specific offsets emphasize that the projected impacts of climate and climate change on near-surface biodiversity and ecosystem functioning are inaccurately assessed when air rather than soil temperature is used, especially in cold environments. The global soil-related bioclimatic variables provided here are an important step forward for any application in ecology and related disciplines. Nevertheless, we highlight the need to fill remaining geographic gaps by collecting more in situ measurements of microclimate conditions to further enhance the spatiotemporal resolution of global soil temperature products for ecological applications